clc;
clear all;
close all;

filename='data\modis\MOD021KM.A2010128.1645.005.2010129013524.hdf';
[image,imageSize] =loadModis1BData(filename);


image=normalize(image);

% modis rgb band number is  [1,4,3]
RGBBand=[1,4,3];
showRGBImage(image,RGBBand);

image=image(1:500,1:500,:);
showRGBImage(image,RGBBand);
label=cluster(image,5,'Kmeans');
figure();
imagesc(label);

